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Recent measurements of the cosmic microwave background (CMB) anisotropies by Planck provide 
a sensitive probe of dark matter annihilation during the cosmic dark ages, and specifically constrain 
the annihilation parameter / e ff (av)/m x . Using new results (Paper II) for the ionization produced 
by particles injected at arbitrary energies, we calculate and provide / e ff values for photons and e + e - 
pairs injected at keV-TeV energies; the / e ff value for any dark matter model can be obtained straight¬ 
forwardly by weighting these results by the spectrum of annihilation products. This result allows the 
sensitive and robust constraints on dark matter annihilation presented by the Planck Collaboration 
to be applied to arbitrary dark matter models with s-wave annihilation. We demonstrate the validity 
of this approach using principal component analysis. As an example, we integrate over the spectrum 
of annihilation products for a range of Standard Model final states to determine the CMB bounds 
on these models as a function of dark matter mass, and demonstrate that the new limits generically 
exclude models proposed to explain the observed high-energy rise in the cosmic ray positron fraction. 

We make our results publicly available at http: //nebel. rc. f as.harvard. edu/epsilon. 

PACS numbers: 95.35.+d,98.80.Es 


I. INTRODUCTION 

Dark matter (DM) annihilation or other new physics 
could inject electromagnetically interacting particles into 
the universe during its early history, with potentially 
wide-ranging observable consequences. In particular, 
injection of ionizing particles during the cosmic dark 
ages will increase the residual ionization fraction, broad¬ 
ening the last scattering surface and modifying the 
anisotropies of the cosmic microwave background (CMB) 
m- Exquisitely sensitive measurements of the CMB by 
a range of experiments (e.g. WMAP [4], SPT [5j [6] and 
ACT |2), and most recently by the Planck satellite [5], 
can thus place robust and model-independent constraints 
on such energy injections (e.g. [9H2Q]). 

For the specific case of DM annihilation, precise con¬ 
straints can be placed under the assumption that the 
power deposited to the gas is directly proportional to that 
injected at the same redshift, with some efficiency factor 
/ e ff. More generally, a model- and redshift-dependent ef¬ 
fective efficiency factor f(z) can be defined as the ratio 
of deposited power to injected power at a given redshift 
(within some arbitrary volume element), even if much of 
the deposited power originates from energy injections at 
earlier times. It has been demonstrated |l6j that given 
a set of f(z ) functions characteristic of Weakly Interact¬ 
ing Massive Particle (WIMP) models for DM [Tl], the 
impact on the CMB is identical at the sub-percent level, 
up to an overall normalization factor, suggesting that ap¬ 
proximating f(z) by a constant / e ff is reasonable. 

In order to convert from deposited power into pertur¬ 
bations to the ionization history, the Planck Collabora¬ 
tion [8] assumed a simple prescription for the fraction 
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of deposited power proceeding to ionization. This pre¬ 
scription was set by the corresponding fraction for 3 keV 
electrons, as determined by detailed Monte Carlo studies 
of the interactions of such electrons by the gas il y 
[24]; we refer to this as the “3 keV” prescription. Ear¬ 
lier works often took an even simpler approach, the so- 
called “SSCK” prescription (based on work by Shull and 
van Steenberg [25], and Chen and Kamionkowski m, in 
which a fraction (1 — x e )/3 of deposited power proceeds 
into ionization; here x e is the ionization fraction (as rele¬ 
vant to the gas species in question; if helium is neglected, 
this is just the hydrogen ionization fraction). In practice, 
these prescriptions do break down; however, since the 
ionizing energy produced by a given energy injection is 
set by the product of the deposition-efficiency curve f(z) 
and the fraction of deposited power proceeding to ioniza¬ 
tion, any errors in the prescription for the latter may be 
absorbed as corrections to f(z). The Planck Collabora¬ 
tion took this approach, employing corrected f(z) curves, 
to set limits for several specific DM models: in this work 
we employ the same constraints, but demonstrate how to 
generalize them to a much broader class of DM models 
(all those where the velocity-weighted annihilation cross 
section ( av ) can be treated as constant during the cosmic 
dark ages). 

In an accompanying article [26], referred to as Paper II, 
we have derived new results for the mapping from energy 
injection into power contributing to ionization, for e + e - 
pairs and photons injected at arbitrary energy (in the 
keV to multi-TeV range) and arbitrary redshift (during 
the cosmic dark ages). Once an energy injection history 
is specified, the power proceeding into ionization at any 
redshift can be determined. In particular, in Paper II we 
have provided corrected f(z) curves for DM annihilation 
(to photons and e + e - pairs) that compensate for the er¬ 
rors in the two most popular ionization prescriptions in 
the literature, “SSCK” and “3 keV”. The only assump- 
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tion on the annihilation model, in deriving these curves, 
is that the redshift dependence of the injected energy is 
dE/dVdt oc (1 + z ) 6 (i.e. the power injected per unit 
volume is proportional to the square of the density), and 
that the spectrum of annihilation products is the same 
at all redshifts. To obtain the f(z ) curve for any given 
model, one need then only integrate over the spectrum 
of photons and e + e - pairs produced by the annihilation. 

In this work we demonstrate that, in accordance with 
previous studies exploring a smaller range of DM models 
mmm, the impact of any conventional DM anni¬ 
hilation model on the CMB can be captured by a sin¬ 
gle normalization factor. The shape of the perturbation 
to the anisotropy spectrum is independent of the spec¬ 
trum of annihilation products, and can be accurately ob¬ 
tained by assuming a constant redshift-independent f(z). 
We derive and present the model-independent weight¬ 
ing function that determines the normalization factor, 
given a (corrected) deposition-efficiency curve f(z). Hav¬ 
ing obtained this weighting function, we can convert the 
corrected f(z ) curves presented in Paper II directly into 
for photons and e + e _ pairs injected by DM an¬ 
nihilation at some arbitrary energy E. Computing f e ^ 
for any conventional model of annihilating DM (i.e. one 
without a novel redshift-dependence) is then simply a 
matter of integrating over the spectrum of annihilation 
products. As an example, we perform this calculation 
for the DM masses and Standard Model final states con¬ 
tained in PPPC4DMID m > and determine the resulting 
constraints from Planck on the annihilation cross sec¬ 
tions for these models. 

We first review, in Sec. [TTJ the interpretation of the 


corrected f(z) curves that we employ. In Sec. Ill 


we 


perform a principal component analysis (PCA) to derive 
the weighting function, construct / e ff as a function of 
injection energy and species, and demonstrate that the 
f e ff results are insensitive to various choices in the analy¬ 
sis. We demonstrate the application of these results to a 
range of DM models in Sec. and then present our con¬ 
clusions in Sec. [V] The appendices provide information 
on supplemental data files, which we make available at 
http://nebel.rc.fas.harvard.edu/epsilon, and re¬ 
view some aspects of PC A and the derivation of the 
weighting function. 


II. THE CORRECTED f(z ) CURVES 

We use the f(z) curves from Paper II [26], and re¬ 
fer the reader to that work for a detailed description of 
their derivation. However, it is worth recapping the def¬ 
inition and properties of the three types of f(z) curves 
we employ, which we will denote f SSCK (z ), f 3keV (z) and 
f sim (z). In the notation of Paper II, these curves cor¬ 
respond to (respectively) f lon i SSCK ( z ^ jion,3keV^ an q 
f sim (z). We will refer to these curves generally as f(z) 
curves. 

The “ion” superscript indicates that the deposition- 


efficiency curves have been corrected to ensure that when 
used with the stated prescription (“SSCK” or “3 keV”), 
one recovers the correct total power into ionization. This 
is the appropriate condition for constraints arising from 
perturbations to the ionization history, which in turn 
modify the anisotropies of the CMB. It has been shown 
that in the case of annihilating DM, other deposition 
channels can essentially be ignored: for example, the lim¬ 
its from high-energy free-streaming photons (discussed 
briefly in e.g. 0 ), distortions to the CMB spectrum 
f28ll3T] and modifications to the gas temperature (e.g. 
PH [32] ) are all much weaker (for s- wave-annihilating 
DM), and ignoring production of Lyman-a photons alto¬ 
gether only changes the constraints by ~ 5% [mm]. ^ c ~ 
cordingly, in this paper we will only ever use f(z) curves 
corrected to obtain an accurate description of the power 
proceeding into ionization (in preference to the other de¬ 
position channels), and the “ion” superscript is unneces¬ 
sary. 


The f sirn (z ) curve is obtained from a simplified esti¬ 
mate for the power deposited into each channel. The 
overall deposited power is corrected by subtracting 
previously-unaccounted power lost into photons below 
10.2 eV, which are no longer sufficiently energetic to 
interact efficiently with the gas (see Paper II [26] for 
details), and then the remainder is multiplied by the 
“SSCK” or “3 keV” prescription to obtain the power into 
ionization. This approach is denoted “approx” in m , in 
contrast to the “best” treatment that yields the deposi¬ 
tion fractions which the f SSCK (z) and f 3keV (z) functions 
are designed to reproduce; it relies on the assumption 
that the “SSCK” or “3 keV” prescriptions are accurate 
once the energy losses of higher-energy electrons (above 3 
keV) to very low-energy photons are taken into account. 
However, as shown for a range of specific models in ng, 
this simplified approach gives results comparable to the 
more detailed calculation. 


To summarize, the / SSCK (z), / 3keV (z) are our best es¬ 
timates for the appropriate corrected f(z) curves to use 
with studies that assumed the “SSCK” or “3 keV” pre¬ 
scriptions, respectively; the f(z) curves have been cor¬ 
rected to take into account errors in those prescriptions. 
The f sim (z) curve is the result of a simpler calculation 
and can be used as a cross-check. We plot all three sets 
of curves in Fig. [lj as taken from Paper II [26]. Note 
that all these curves only account for the energy injected 
by annihilation of the smooth component of DM with 
the cosmological average density; they do not take into 
account structure formation at low redshifts. However, 
it has been shown that the contribution from structure 
formation to the CMB constraints is generally very sub¬ 
dominant for 5-wave annihilation m- 
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FIG. 1: Corrected f(z) functions for particles injected by DM annihilation, as a function of injection energy and redshift of 
absorption. In the left panel we use the “3 keV” baseline ionization fractions (so these f 3keV (z) curves should be used with 
analyses that employed the same prescription); in the center panel we use the “SSCK” baseline. In the right panel we plot the 
simplified channel-independent f sim (z) curve. The upper row describes e + e _ pairs (the x-axis “energy” label here indicates 
the kinetic energy of a single member of the pair at injection), the lower row describes photons. 


III. CONVERTING CORRECTED f(z) CURVES 
TO CMB BOUNDS 

The most recent and strongest constraints on DM an¬ 
nihilation from CMB anisotropies come from an analysis 
of Planck data [8] . These constraints were computed us¬ 
ing the “3 keV” baseline prescription, and were expressed 
in terms of a bound on f e s(crv) /m^M- For the models 
tested by those authors, / e ff was determined by evaluat¬ 
ing the appropriate corrected f(z ) curve for the model 
at z = 600, as originally suggested and validated in Em 
Analogous f e ^ values for arbitrary photon and e + e en¬ 
ergies can be read off directly from Fig. [l] using the 
“3 keV” baseline prescription or the simplified f sim (z ) 
curves (as discussed above, the latter approach is simpli¬ 
fied and likely to be less accurate, but it gives a rough 
estimate of the possible systematic errors arising from 
details of the first method). 

A more careful derivation of / e ff - that is, the constant 
value of f(z) which would have the same impact on the 
CMB as the true corrected f(z) curve - can be achieved 
by PC A, following the approach of m - 1 We review the 


1 For convenience, in this work we use the marginalized Fisher ma¬ 
trices derived and made publicly available by m-, while those 
results employed WMAP measurements of the cosmological pa- 


essentials of the Fisher matrix and PC A in Appendix [Bj 
Performing PCA on a range of f(z) curves correspond¬ 
ing to assorted DM models, HH [2D] found that the vast 
majority of the variance - exceeding 99.9% - was con¬ 
tained in the first principal component, peaking around 
z ~ 600. Related studies on Hang independently de¬ 
termined, using different methods, that the signal was 
largely controlled by a single parameter, set by the be¬ 
havior of f(z) at z ~ 600. 

As described in [j2(|, the “dot product” of a given f(z) 
curve with the first principal component thus completely 
determines its detectability. This dot product can be 
re-expressed as an integral over dln(l + z), with the inte¬ 
grand being the product of the f(z) curve with a weight¬ 
ing function W(z). The weighting function can be deter¬ 
mined from the first principal component, as described 
in (16j 20 ; we provide a detailed derivation in Appendix 
[Cl (An alternate PCA, based on the ionization history 
rather than the energy deposition history, was performed 
in [34].) 

In principle, the weighting function W(z) depends on 
the details of the experiment (specifically, its sensitivity 


rameters to determine the fiducial cosmological model, we expect 
the impact of updating to Planck cosmological parameters to be 
extremely small. 
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to different multipoles £, in both temperature and polar¬ 
ization). It will also depend on the ionization prescrip¬ 
tion, as this affects the mapping between f(z) and the 
CMB. In practice, however, the impact of these choices 
turns out to be very small. 

We performed a new PCA, building our model space 
from the f(z) curves for annihilation to e + e _ and pho¬ 
tons, for 20 log-spaced energies between 1 keV and 10 
TeV (for 40 models total). This is considerably more 
general than previous studies which focused on the GeV- 
TeV WIMP parameter space, but again we find that the 
first principal component overwhelmingly dominates the 
variance (accounting for more than 99.7% in all cases). In 
other words, the space spanned by the possible imprints 
of conventional 8-wave DM annihilation on the CMB is 
(within the confines of the approximate Fisher-matrix 
approach) essentially one-dimensional. 

In order to justify computing the constraints on con¬ 
stant f(z ) and applying them to arbitrary annihila¬ 
tion models, we need to show that scenarios with con¬ 
stant f[z ) also produce a perturbation to the CMB 
anisotropies lying within this one-dimensional space. To 
test this question, we consider two energy injection his¬ 
tories, one with f(z ) = constant, the other with f(z) 
given by the first principal component, with the same 
DM mass and annihilation cross section. Suppose the 
f(z) curves are normalized such that they correspond to 
the same signal significance. As discussed in Appendix 
|B|(and using the notation defined there), if v\ and v 2 are 
the discretized version of the f(z) curves, this condition 
implies vf F z v\ = F z n 2 , where F z is the Fisher matrix 
(as determined in nsp- 

The significance of the difference between these energy 
injection histories, normalized to the overall significance 
of either history, can be estimated within the Fisher ma¬ 
trix approach as: 




l(v 2 - vi) t F z (v 2 — Vi) 


v[F 2 v 1 



( 1 ) 

Here we have used the normalization condition above, 
and also the symmetry of the Fisher matrix to write 
v[F z V2 = V 2 F z v\. We computed the bracketed quan¬ 
tity on the RHS of Eq. [I] taking fi(z) to be constant 
and f 2 (z) to be the first principal component. We found 
that in all cases the bracketed quantity was below 0.001, 
indicating that to a good approximation, the first princi¬ 
pal component also captures the effect of constant f(z). 
Thus it is reasonable to derive constraints on the case of 
constant f(z) and convert them to constraints on general 
DM annihilation histories using the weighting function. 

We repeated the analysis for each independent combi¬ 
nation of the following choices: 


• A Planck-\ike experiment vs a cosmic variance lim¬ 
ited (CYL) experiment (as defined in p"6]). 


• Including multipoles up to £ = 2500, vs i = 6000. 


• Using f SSCK (z ), f 3keV (z) and f sirn (z ) to generate 
the space of f(z) curves. 

• Choosing “SSCK” vs “3 keV” as our baseline ion¬ 
ization prescription (to translate the f(z) curve 
into its impact on the CMB) when using the f sim (z ) 
curve. 

In all cases we used a Fisher analysis to estimate de¬ 
tectability in the CMB, as described in m- We normal¬ 
ized the weighting function as described in [20] and de¬ 
rived in detail in AppendixlQ so that f W(z)d\n(l-\-z) = 
1. Thus for any arbitrary f(z), f f(z)W(z)dln(l + z) 
measures detectability relative to the case of f(z) = 1. 

There are visible differences in the weighting func¬ 
tions thus derived, with the largest impacts coming from 
changing the sensitivity to different multipoles - that is, 
by considering a Planck-like experiment, vs a CVL ex¬ 
periment up to i — 2500, vs a CVL experiment up to 
i = 6000. However, the resulting variation in f e ^ com¬ 
puted by the different W (z) functions is always negligible 
(percent-level or less). We show the various W(z) curves 
in Fig. [2j in Fig. [3] we display the impact on the f e & 
fractions of scanning over all derived W(z) curves. 

We recall that if the f SSCK (z) and f 3keV (z) curves 
are each used with their specified ionization prescription, 
they correspond to identical changes to the ionization his¬ 
tory. Consequently, their impacts on the CMB must be 
the same (if contributions from non-ionization channels 
are ignored). However, since the f(z) curves are differ¬ 
ent, one might expect that the weighting function W(z) 
will be different - and so must be chosen to match the 
ionization prescription under which it was derived - in 
order to compensate. We discuss this point in detail in 
Appendix[Dj but it turns out that in practice, the similar¬ 
ity of the different ionization prescriptions and the unit 
normalization of W(z) means the effect is tiny. Instead, 
the difference in baseline ionization prescriptions is cap¬ 
tured in the / e ff parameter, and the derived constraints 
on that parameter; for example, a higher baseline pre¬ 
scription for ionization would correspond to lower f(z) 
for all models, and hence lower / e ff, but the constraints 
on a given / e fr would be stronger (since the higher ion¬ 
ization prescription corresponds to a larger effect on the 
CMB). 

We have explicitly tested the effect of using the 
“wrong” weighting function - i.e. using W(z) derived 
using the f 3keV (z) curves, in order to estimate / e ff for 
/ ssck ^) curves> The effect is included in the bands in 
Fig-! and is negligible. 

Consequently, for the purposes of computing / e ff, it is 
entirely adequate to use a single weighting function; we 
choose the weighting function derived for a CVL exper¬ 
iment including t up to 2500, built from the f 3keV (z) 
curves. We note that this choice is largely arbitrary, as 
the resulting value of / e ff is very stable. This weight¬ 
ing function is likely to be broadly applicable beyond the 
case of DM annihilation, although it may not capture 
all the variance for energy injection histories that dif- 
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FIG. 2: Weighting functions for the f(z) curves, derived by PCA. Solid lines show the results for a CVL experiment covering 
£ = 2 — 2500, dashed lines show the results for a CVL experiment covering £ = 2 — 6000, dotted lines show the result for a 
Planck -like experiment covering £ — 2 — 2500. Black lines use the first procedure for estimating ionization history, red lines use 
the second procedure with “3 keV” baseline prescription, blue lines use the second procedure with “SSCK” baseline prescription. 
Left panel: suitable for “3 keV” baseline. Center panel: suitable for “SSCK” baseline. Right panel: “universal” weighting 
function recommended for general use. 


fer markedly from that of conventional DM annihilation. 
We leave a more detailed PCA, based on a larger range 
of possible injection histories, for future work. 

One can integrate any electron and photon spectra 
produced by DM annihilation over the f e ff curves pre¬ 
sented in Fig. [3j in order to determine / e ff for an ar¬ 
bitrary model. We show results separately for the “3 
keV” and “SSCK” baseline prescriptions; however, when 
making comparisons with the Planck constraints [8], the 
“3 keV” baseline prescription should always be used. We 
also show the effect of taking our second simplified proce¬ 
dure for computing / corr (z), and the effect of evaluating 
/corr(600) rather than using the weighting function; both 
effects are rather small, at the < 10% level. 

The /efr curves (Fig. [3| exhibit significant structure as 
a function of energy. At high energies, both e + e - pairs 
and photons asymptote to / e ff « 0.4; this convergent be¬ 
havior is expected, since at high energies electrons and 
photons both participate in a pair production / ICS cas¬ 
cade. Both species demonstrate some minor structure 
in the GeV-TeV range but generally maintain / e ff ^ 0.4, 
with a peak of / e ff ~ 0.7 for electrons at injection energies 
around 10 GeV. Both species then exhibit a large bump 
in / e ff, peaking at energies around 10-100 MeV, and a 
trough at O(MeV) energies. This trough occurs because 
mildly relativistic electrons suffer large losses into con¬ 
tinuum photons with energies below 10.2 eV (see Paper 
II [26] for an in-depth discussion). Accordingly, MeV- 
scale electrons have small / e ff, and the dominant energy 
loss process for MeV-scale photons is Compton scatter¬ 
ing, which in turn produces more mildly-relativistic elec¬ 
trons. The peak at slightly higher energy occurs because 
10-100 MeV electrons lose most of their energy by inverse 
Compton scattering into 0(100 —10 4 ) eV photons, which 
are efficient ionizers. 

At low energies, below ^ 1 — 100 keV (depending on the 
redshift), photons become efficient photoionizers, causing 
a large / e ff. In contrast, the e + e~ / e ff at low kinetic ener¬ 
gies is rather low, of order 0.3 - but as noted in m, this is 


because in the context of DM annihilation, almost all the 
injected power is bound up in the mass energy of the pair 
in this case. Consequently, the absorption of the kinetic 
energy is almost irrelevant: what matters is the absorp¬ 
tion of the photons produced when the positron annihi¬ 
lates. From inspection of / e ff for ~ 0.5 MeV photons, 
one would expect / e ff of 0.3 for pairs produced nearly at 
rest, consistent with the results of the calculation. 

The primary systematic uncertainties in this calcula¬ 
tion arise from the choice of f(z) curves (/ SSCK >3keV^ or 
/ sim (z)). Both of these curves involve an approximation; 
on one hand matching the power into ionization exactly 
(up to the limitations of our approximate treatment of 
the low-energy photons) while not attempting to match 
the heating and Lyman-a channels, on the other approx¬ 
imating all channels by the cooling of 3 keV electrons 
after accounting for energy lost to continuum photons. 
Some approximation of this form is inevitable, given that 
the true relative losses into the different channels are not 
independent of the injected particle energy as tacitly as¬ 
sumed by both the “SSCK” and “3 keV” prescriptions. 
As discussed more extensively in [18], we estimate the 
potential systematic error due to these approximations 
and our simplified treatment of the low-energy photons 
to be at the ^ 10% level. For some channels, the con¬ 
straints we present may also be slightly too conservative, 
i.e. weaker than the true limits, at a similar ^ 10% level, 
due to our neglect of energy deposition by protons and 
antiprotons [55] . 


IV. CONSTRAINTS FOR SIMPLE 
ANNIHILATION CHANNELS 

In this section we apply the results of this work 
to the Standard Model final states provided in the 
PPPC4DMID package [27], following the earlier analysis by 
[315]. PPPC4DMID provides the fluxes of positrons and pho¬ 
tons produced by DM annihilation to 28 Standard Model 
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FIG. 3: / e ff coefficients as a function of energy for e + e _ (left column) and photons (right column); appropriate for analyses 
with baseline prescription “3 keV” (top row) or “SSCK” (bottom row). The widths of the red and green bands indicate the 
impact of scanning over all the derived weighting functions (see Fig. f2|) . Red stars are derived from the f SSCK (z) or f 3keV (z) 
curves as appropriate, green stars from the f sim (z) curves. Diamonds indicate / e ff evaluated as /(600), rather than using a 
weighting function, and are labeled “simple photon-loss rescaling” in the legend. 


final states, for masses ranging from 5 GeV to 100 TeV. 
The fluxes include electroweak final state radiation m , 
which can be crucial at high DM masses; thus, for exam¬ 
ple, there is a small but non-zero f e Q for annihilation to 
z/zq arising from the electroweak corrections. 

We only show results up to 10 TeV, as above this 
mass, additional processes not considered in m may 
contribute to the cooling of electrons and photons in the 
early universe. However, we expect that the f e ff results 
will be similar at higher DM masses, since the processes 
neglected in HQ will only modify the details of the pair 
production cascade experienced by particles injected at 
high energy, which occurs much faster than a Hubble time 
and so is anticipated to have little effect on the fraction 
of power eventually deposited. 

We take the positron and photon spectra provided in 
PPPC4DMID, 2 assume that by charge symmetry the spec¬ 
trum of e~ is equal to that of e + , and integrate the re¬ 
sulting e + e - and photon spectra over the f e s(E) curves 
presented in Fig. [3j Specifically, if a given annihilation 
produces a spectrum ( dN/dE) e + of positrons, and a spec¬ 
trum (dN/dE)ry of photons, the weighted f e ff(m x ) curve 


2 http://www.marcocirelli.net/PPPC4DMID.html 


is given by: 


fes(m x ) = 


Jo”* EdE 2/; 


/ e 1 e 
eff 


'(E) (45H 


dE ) e+ 


f: s (E)d it 


2 


( 2 ) 

where E is the energy of the positron or photon (the 
factor of 2 comes from the electron spectrum). 

We use the best-estimate f e fi(E) curves suited for 
the “3 keV” baseline prescription, as we wish to ap¬ 
ply the constraints derived by the Planck Collabora¬ 
tion [ 8 ]. We take the Planck limit to be f e ^(av)/m x < 
4.1 x 10 -28 cm 3 /s/GeV, which is the bound set by the 
Planck temperature and polarization data. The limit 
can be strengthened by ^ 15% by the inclusion of either 
lensing or external datasets. We show both the f e ^(m x ) 
curves for each channel and the resulting constraints on 
the annihilation cross section in Fig. [4] For DM masses 
below 5 GeV, we show estimated limits for DM anni¬ 
hilation to e + e - , 77 or VV —>■ 4e - the last channel 
corresponds to DM annihilating to a pair of intermedi¬ 
ate vector bosons, each of which subsequently decays to 
e + e - . The results for the first two channels are taken 
simply from Fig. [ 3 ) whereas for VV 4e, the e + e - 
spectrum is taken to be constant in dN/dE between the 
kinematic limits m e < E < m x — m e . For both channels 
producing e + e - , there will also be an FSR contribution, 
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FIG. 4: The upper panel shows the / e ff coefficients as a function of DM mass for each of a range of SM final states, as indicated 
in the legend. The VV —>■ 4X states correspond to DM annihilating to a pair of new neutral vector bosons V, which each 
subsequently decay into e + e _ , or t + t~ (labeled by X ). The lower panels show the resulting estimated constraints from 

recent Planck results [ 8 ], as a function of DM mass, for each of the channels. The left panel covers the range from keV-scale 
masses up to 5 GeV, and only contains results for the e + e _ , 77 and VV —> 4e channels; the right panel covers the range 
from 5 GeV up to 10 TeV, and covers all channels provided in the PPPC4DMID package m- The light and dark gray regions 
in the lower right panel correspond to the 5 a and 3 cr regions in which the observed positron fraction can be explained by DM 
annihilation to for a cored DM density profile (necessary to evade 7 -ray constraints), taken from [5BJ. The solid yellow 

line corresponds to the preferred cross section for the best fit 4-lepton final states identified by [ 33 , who argued that models 
in this category can still explain the positron fraction without conflicts with non-observation in other channels. The red and 
black circles correspond to models with 4e (red) and 4 /i (black) final states, fitted to the positron fraction in [38] : as in that 
work, filled and open circles correspond to different cosmic-ray propagation models. 


but its effect is generally small (at the percent level). 

In general, we see that the final states considered fall 
into three categories: 

• Final states where the bulk of the power pro¬ 
ceeds into e + e - and photons, where at masses 
above 100 GeV the constraint approaches (av) < 
10 _27 (ra x /lGeV) cm 3 /s. 

• Annihilation to neutrinos, where the constraint 
arises entirely from electroweak corrections, and is 
negligible below ^ 200 GeV; at O(TeV) masses, 
cross sections as low as a few xlO -23 cm 3 /s can be 
constrained. Interestingly, this bound is competi¬ 
tive with that placed by IceCube from observations 
of galaxy clusters S 3 , the Galactic Center [42] , and 
the Milky Way halo [43], and unlike those limits is 
independent of uncertainties in the local DM den¬ 
sity, the DM distribution, and the amount of DM 


substructure. 

• A band with a width of roughly a factor of 150% in 
(av) that encompasses all the other channels stud¬ 
ied, which at high masses corresponds to (av) < 
2 — 3 x 10 _27 (ra x /lGeV) cm 3 /s. 

Accordingly, for any linear combination of these final 
states that does not contain a significant branching ratio 
for DM annihilation directly to neutrinos, one must have 
(av) < 3x 10 _27 (ra x /lGeV) cm 3 /s. It is thus challenging 
to obtain the correct thermal relic cross section for s-wave 
annihilating DM with mass much below m x ~ 10 GeV, 
without violating these limits (although models with sup¬ 
pressed annihilation at late times may still be viable, 
e.g. asymmetric DM models or the scenarios proposed in 
[441145] ). At higher masses, the cross sections constrained 
are well above the thermal relic value, but are highly rele¬ 
vant for DM explanations of the positron excess observed 









































by PAMELA [46], Fermi [47] and AMS-02 [48]. For ex¬ 
ample, m identified a favored scenario by which to ex¬ 
plain this excess via DM annihilation, respecting all ex¬ 
isting constraints: 0.5-1 TeV DM annihilating through a 
light mediator to a mixture of 75% t + t - and 25% e + e~. 
The required cross section is roughly (av) ~ 5 x 10 -24 
cm 3 /s at 500 GeV and roughly ~ 1.5 x 10 -23 cm 3 /s at 1 
TeV. This is clearly above the bounds for the contribut¬ 
ing annihilation channels in Fig. [4] Similarly, [38] found 
that only annihilation through a light mediator to muons 
or electrons could explain the excess without violating 
bounds from studies of dwarf galaxies in gamma rays. 
The preferred cross section was ^ 1.3 —1.6 x 10 -24 cm 3 /s 
at a DM mass of ~ 350 GeV for the 4e final state, and 
~ 6 — 8 x 10 -24 cm 3 /s at a DM mass of ~ 600 — 750 GeV 
for the 4/i final state. Again, these preferred points are 
well above the bounds we derive for the corresponding 
channels. 

Fig. El shows (1) the cross section band favored for 
the models of as a solid yellow line, (2) the points 
favored by [38j as circles in red (4e final state) or black 
(4/i final state), and (3) a gray contour describing the 
region in which DM annihilation to fi + fi~ can explain 
the AMS-02 data [36]. These favored regions should be 
compared, respectively, to (1) an appropriately weighted 
combination of the green dot-dashed and red dot-dashed 
lines in Fig. 4] (it is a conservative choice to take the 
green dot-dashed line to set the limit), (2) the red dot- 
dashed and black dot-dashed lines, and (3) the solid black 
line. We see that in general these models are ruled out by 
more than a factor of 2, and up to an order of magnitude. 
The models presented in [45 were studied in [8], and are 
similarly ruled out by these constraints. 

V. CONCLUSIONS 

In this work, we have mapped out the effec¬ 
tive deposition efficiency, / e ff, for photons and e + e - 
pairs injected into the early universe by DM an¬ 
nihilation, at O(keV-TeV) kinetic energies. We 
have made our key numerical results available at 
http://nebel.rc.fas.harvard.edu/epsilon. 

We have demonstrated that it is sufficient to use a sim¬ 
ple universal weighting function to convert any redshift- 
dependent energy deposition history (originating from 
conventional s- wave-dominated DM annihilation) into a 


redshift-independent efficiency factor, confirming earlier 
results for the WIMP regime. This efficiency factor then 
converts Planck limits on energy injection from the CMB 
into constraints on the DM model in question. Any DM 
model featuring dominantly s-wave annihilation, with 
the only redshift dependence in the annihilation rate 
arising from the DM density squared, can thus be con¬ 
strained using these data. We emphasize that this does 
not subsume the limits presented by the Planck Collab¬ 
oration, but confirms the validity of the procedure em¬ 
ployed in that analysis, and demonstrates how to extend 
those bounds to more general DM models. 

We have performed an example calculation for 28 dif¬ 
ferent annihilation channels over the DM mass range 5 
GeV - 10 TeV, and demonstrated that except for an¬ 
nihilation directly into neutrinos, at high masses the 
constraint from the CMB generically lies in the range 
(av) < 1 — 3 x 10 _27 (ra x /lGeV) cm 3 /s. Even for anni¬ 
hilation to neutrinos, at TeV-scale masses we can set an 
upper bound on the cross section of order a few times 
10 -23 cm 3 /s, comparable to constraints from IceCube 
observations of the Milky Way and galaxy clusters; this 
bound arises from electroweak final state radiation pro¬ 
ducing e + e - pairs and photons. 

The constraints we obtain appear to rule out annihilat¬ 
ing DM models proposed to explain the observed excess 
of cosmic-ray positrons, including models which could 
not be excluded by any previous studies. This constraint 
could be evaded if the DM annihilation is suppressed at 
low velocities or early times. 
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Appendix A: Supplementary Materials 

We make available 3 a selection of .fits and .dat files containing the f e ^ values presented in Figs. [3] and [I] (with 
the exception of the VV —>■ 4e channel below 5 GeV, which can be trivially obtained - at least at the level we have 
calculated it - from the results for e + e - ), as well as the universal annihilation weighting function shown in Fig. [ 2 J 


3 http://nebel.rc.fas.harvard.edu/epsilon 




9 


We also provide an example Mathematica notebook to demonstrate the use of the .fits files. 

The .fits files are as follows: 

feff_summary_species_base.fits 

Here “species” can be “elec” (representing e + e~ pairs) or “phot” (representing photons), and “base” can be “SSCK” 
or “3 keV”, corresponding to the choice of baseline ionization prescription. For example, the “3 keV” files should be 
used with the constraints on DM annihilation of [8]. Each file contains the following arrays: 

• REDSHIFT: this 63-element array provides the abscissa for deposition redshift. 

• WEIGHTFN: this 63-element array provides the weighting function shown in the third panel of Fig. [2j sampled 
at the redshifts given by the REDSHIFT array, and normalized so that fW(z)d ln(l + z) = 1 (up to numerical 
error). 

• L0G10ENERGY: this 40-element array gives log 10 (E/eV), where E is the kinetic energy of one member of the pair 
of injected particles (for e + e - pairs) or the energy of a single injected photon. 

• FEFF_BEST: this 40-element array gives / e ff computed using the appropriate f hase (z) curve, weighted by the 
weighting function. 

• FEFF_RESCALING: this 40-element array gives / e ff computed using the appropriate f sirn (z ) curve, weighted by 
the weighting function. 

• FEFF_BEST_600: this 40-element array gives f e ff computed using the appropriate f hase (z) curve, evaluated at 
z = 600. 


feff_PPPC_3keV.fits 

This file contains the f e ff curves for different DM masses and annihilation channels as shown in Fig. [4j at DM 
masses above 5 GeV. Below 5 GeV, the results for the e + e - , 77 and VV 4e channels are easily obtained from the 
feff_summary_species_3keV.f its files. The arrays in this file are: 

• DM_MASS: this 62-element array gives the DM mass in GeV. 

• CHANNEL: this 28-element array labels the channels as shown in the legend of Fig. |4j 

• FEFF: this 28 x 62 array contains the / e ff values for each DM mass and channel, derived using the appropriate 
f 3keV (z) curves and the weighting function. 

The .dat files simply contain condensed summaries of the .fits files, and are as follows: 


feff_PPPC_3keV.dat 

As feff_PPPC_3keV.fits, except that the / e ff values are laid out in a grid, with columns corresponding to different 
channels (first row describes the different final states) and rows corresponding to different DM masses (first column 
gives the DM mass in GeV). 

f ef f _base . dat 

Contains a subset of the information in feff .summary _electron_base. fits and feff .summary _photon_base . fits: 
the first column is the L0G10ENERGY array from these files, describing the injection energy; the second is the best- 
estimate /eff value (computed using the f hase (z) curves and the weighting function) for e + e - pairs; the third is the 
best-estimate / e ff for photons, “base” can take the values “SSCK” or “3 keV” and should be chosen based on the 
ionization prescription used in the analysis to which comparison is desired. 

we ighting_f unction. dat 

Contains the weighting function W(z) shown in the third panel of Fig. [ 2 J sampled at the redshifts 1 + z given in 
the first column. 
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Appendix B: Review of Principal Component Analysis 

Suppose we have some arbitrary energy deposition history, characterized by a (corrected) efficiency function f(z). 
Let us discretize the f(z ) function into n z bins centered at Zi, i = 1 ..n z , thus describing the energy deposition by a 
column vector v = (/(^i), f[z^), f(z 3 ), f(z nz )). We follow [16] in treating the energy deposition profile f(z ) in each 
bin as a Gaussian in ln(l + z), normalized such that its integral with respect to dln(l + z) is equal to the logarithmic 
bin width A ln(l + z). 

Within the approximation of linearity (which is a good approximation for energy injections at the level constrained 
by Planck [16]), the impact of this deposition history on the anisotropies of the cosmic microwave background can be 
written in the form, 


h = A C t = Tv, 


(Bl) 


where T is some transfer matrix of dimension ni x n z , and n£ is the number of multipoles we consider. Following 
[16], in our notation each element of T (Ta) is itself a three-element vector, holding the perturbations to the TT, TE 
and EE anisotropy spectra at that £. T can be computed element-by-element, as described in m, by taking Vj such 
that Vi = Sij, and computing the effects on the CMB using RECFAST [50 (or a similar code such as CosmoRec m or 
HyRec [52] ) and CAMB [53]. Similarly, we can determine a corresponding x 6 transfer matrix describing the effect of 
varying the standard six cosmological parameters. 

Once T has been mapped out, we can construct the n z x n z Fisher matrix F e as, 


{T e )ij — Tgi • E^ Tij, 


(B2) 


where E^ is the appropriate covariance matrix for the anisotropy spectra (see m for the explicit form). Schematically, 
F e = T t E - 1 T; it is a symmetric matrix, due to the symmetry of E, and its diagonal elements ( F e )u describe the 
(squared) signal significance per energy deposition at redshift Zi, before marginalization over the existing cosmological 
parameters. 

To perform this marginalization, one first constructs an expanded Fisher matrix Fq, using Eq. B2 but expanding 
T to include the effect of the standard six cosmological parameters. We can write: 


Fn = 


F e F v 
FT F r 


(B3) 


where F e is the pre-marginalization Fisher matrix as defined above, F c is the Fisher matrix for the cosmological 
parameters only, and F v describes the cross terms. Then the marginalized Fisher matrix, describing the detectability 
of an energy injection after projecting out the perturbations parallel to the standard cosmological parameters, is given 
by F z —F e — FyF^F^ (see [16] for a discussion). 

The principal components are given by the eigenvectors of F z , which we will denote {Pi}, i = l..n z . By construction, 
they form an orthogonal basis for arbitrary f(z) histories, and we can choose them to be ortho normal; as such, one 
can write the vector corresponding to any arbitrary f(z) in the form v = ^2aiPi, where Pi is the ith principal 
component and oti = v • Pi. Each principal component perturbs the CMB anisotropies with significance proportional 
to the square root of the corresponding eigenvalue (to the degree that the approximations inherent in the Fisher 
matrix analysis are valid). Consequently, if the principal components are ranked in order of their eigenvalues, the 
detectability of an arbitrary energy deposition history can be described approximately by aq = v • Pi, the dot product 
with the first principal component. In this sense, Pi acts as a redshift-dependent weighting function that determines 
the detectability of an arbitrary energy deposition history. The accuracy of this “weighting function” approximation 
will be controlled by the fraction of the total variance described by the first eigenvalue (it may also break down for 
specific energy deposition histories where v is accidentally nearly orthogonal to Pi). 


Appendix C: Deriving the Universal Weighting Function for Annihilating DM 

Rather than working in terms of energy deposition histories localized in redshift, one can perform a similar PC A 
starting with any sample of f(z) functions (see [16] for a more complete discussion); the resulting principal components 
will be applicable to any energy deposition history lying within the space spanned by these basis functions. For the 
purposes of setting CMB constraints, any conventional DM annihilation model can be described as a linear combination 
of photons and e + e _ pairs injected over a range of energies. 
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Consequently, we can perform a PC A applicable to all conventional DM annihilation models, by taking as our 
basis functions the f(z ) curves plotted in Fig. [l] We choose 20 log-spaced injection energies for photons and e + e - 
pairs, thus giving rise to 40 f(z ) curves as our sample. As previously, we can sample these curves at n z redshifts, 
translating them into n^-element column vectors, which we will denote U{ for i = 1..40. (More generally, we will 
denote the number of models included in the PC A as tim • in this case um =40.) Let U be the n z x um matrix 
with the ith column given by iq. A general DM annihilation model can be expressed as r = Ua , where a now 
describes the coefficients of the ( not orthogonal) basis vectors iq. Thus we can describe the matrix U as mapping 
from ?tm- dimensional “coefficient space” into n^-dimensional “redshift space”. 

Following the same procedure as previously, we can compute the impact on the CMB of each of the Ui histories, 
thus constructing a ri£ x um transfer matrix Tm that maps from “coefficient space” to the anisotropies of the CMB. 
In terms of the earlier transfer matrix T, Tm = TU. Let us label the associated um x um Fisher matrix (after 
marginalization over the standard cosmological parameters) as Fm, and the expanded pre-marginalization Fisher 
matrix as F 0 m- 

Since the transfer matrix for the cosmological parameters is unchanged, writing Tm in terms of T yields the following 
relation between Fq and Fqm : 


Fqm 


U T F e U U t F v \ 
(U t F v ) t f c ) • 


Consequently, the marginalized Fisher matrix becomes, 


F m = U T F e U - U t F v F~ 1 F^U = U t F z U. 


(Cl) 


(C2) 


Diagonalizing Fm now yields a set of principal components specialized to the case of DM annihilation; these principal 
components are nM-element vectors, and describe coefficient vectors a for the original basis f(z) curves. Accordingly, 
they are orthonormal in “coefficient space”, not “redshift space”. If we denote these principal components as { pi }, 
i = 1then following the reasoning above, we can write a = JT(a ' Pi)Pi ? an d the detectability of the energy 
deposition history characterized by the coefficient vector a is determined approximately by a • pi . 

Due to the similarity of the f(z) curves for different energies and different species, when the eigenvalues are 
calculated, the first principal component accounts for more than 99% of the variance. Accordingly, having specialized 
to the case of DM annihilation, the “weighting function” approximation is expected to be accurate at the sub-percent 
level. 

However, we would prefer to express this dot product directly in terms of the energy deposition history f(z ), or 
its discretized counterpart F, rather than its coefficients with respect to a somewhat arbitrary basis. To achieve this, 
note that since v = Ua, we can write, 


U t F z v = U T F Z Ua = F M ^>2(d-Pi)Pi = ^ A*(a ■ pl)pl 


(C3) 


2=1 2=1 

Here A i is the ith eigenvalue of Fm, and so FmV% — \v% by definition. Taking the dot product with pi then yields, 

p 1 -(U T F z v) = X 1 (a-p{), (C4) 


and consequently, 


a-Pi 



^ (Ufi) T F z ^ 


• V. 


(C5) 


If a constraint on some deposition history fo(z) (lying within the space spanned by conventional DM annihilation) 
has been computed, therefore, the constraint on any other f(z) (at the same DM mass and cross section) can be 
computed by rescaling the bound by a factor: 


(^)-O 


(Up^F, 

((UPi) T F z )-Zo 


(C6) 


A higher f e ^ corresponds to a more detectable f(z ), and hence to a stronger bound. The quantity in square brackets 
now acts as a weighting vector W which converts v to f e ^. In principle, we could use any vector vq to determine the 
reference bound, but in practice the usual choice is f(z) =■1, so / e ff represents “effective deposition efficiency” relative 
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to the case corresponding to total prompt energy deposition (with one of the simplified ionization prescriptions). This 
choice yields: 


W = 


0 Upi) T F z 


(C7) 


We can define a continuous weighting function W ( z ) by interpolation, fixing: 


W(zi) 


(Wh 

d ln(l + z)' 


(C8) 


Then the normalization condition ^2 i (W)i = 1 translates into JT W(zi)d\n(l + z) = 1, and so when the continuous 
limit is taken, f W(z)d\n(l + z) = 1. We can likewise rewrite f e g in terms of the continuous functions f(z ) and W(z ), 


/eff = E /(^i)(W r )< = E f( z i)W(zi)d ln(l + z) 

i i 

J /(z)W(z)dln(l + z). 


This is the weighting function provided in the main text. 


(C9) 


Appendix D: Expected Change to the Weighting Function Due to Ionization Prescription 

As discussed in Sec. m the weighting function has some dependence on the presumed baseline ionization pre¬ 
scription, even though the f SSCK (z) and f 3 keV (z) curves correspond to identical physical ionization histories (when 
properly combined with the appropriate ionization prescriptions). There is also a separate dependence on the pre¬ 
sumed baseline ionization prescription when the f sim (z ) curve is used, but in this case the dependence is physical; it 
affects the mapping from DM models into the CMB, by changing the transfer matrix Tm (and T). 

In the notation of Appendix |c| when f SSCK (z) or f 3 keV (z) curves are used, the Fm Fisher matrix (and resulting 
principal components) are unaffected by the choice of ionization prescription, as these depend only on the (physical) 
CMB signatures of particular DM annihilation models. However, the U and F z matrices, which involve mappings 
between DM models and f(z) curves, or f(z) curves and the CMB, are affected. 

Let U\ and U2 denote the matrices that hold the f(z) curves for two different baseline ionization prescriptions, and 
let 22 ( 2 ) describe the fraction of deposited power proceeding into ionization (as a function of redshift) in those 

two prescriptions. Let I\ and I2 define the n z x n z diagonal matrices such that (. I a )jk = ia( z j) 3 jk- Then since by the 
definition of f hase (z) curves, the power into ionization at any redshift is invariant, we can write I\U\ = I2U2 • 

Now if F z 1 and F Z 2 are the Fisher matrices corresponding to the two baseline ionization prescriptions, we can write 
F m = UjF zl Uy = UlhI^ l F zl I- 1 I l U l = Ujl 2 IpF zl IpI 2 U 2 = U?F z 2 U 2 . It follows that, 

F z2 = I 2 IpF zl IpI 2 . (Dl) 

Then we can immediately write down the relationship between the weighting vectors corresponding to the two pre¬ 
scriptions, 


W 2 


(U 2 Pi) T F z2 

Ei mpi ) T F Z 2) i 

{IpIilhpi) T l2lpF zl IpI 2 

E,: ((IphUypi) T hIpF zA I-'I 2 ) i 


{u^Yf, 


17 l h 


Z 1 J 1 


E« {(UlPlVF; 


IF 1 19. 


z l 1 ! 




T F Z 1) 


'Ei ((UmVF^ii 


-1 



(D2) 


We see that if I\ and I2 only vary by a redshift-independent normalization factor, the effect on W cancels out 
entirely, and this is a quite a good approximation at all redshifts of interest for the signal (see e.g. Paper II [26]). 
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Thus while the redshift-dependence of I 1 does slightly distort the shape of the weighting function, the impact on 
the extracted / e ff is negligible (as demonstrated explicitly in Sec. 


III). 
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